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I/Q Imbalance Aware Widely-Linear Receiver 
for Uplink Multi-Cell Massive MIMO Systems 
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Abstract 

In-phase/quadrature-phase (I/Q) imbalance is one of the most important hardware impairments 
in communication systems. It arises in the analogue parts of direct conversion radio frequency (RF) 
transceivers and can cause severe performance losses. In this paper, I/Q imbalance (IQI) aware widely- 
linear (WL) channel estimation and data detection schemes for uplink multi-cell massive multiple-input 
multiple-output (MIMO) systems are proposed. The resulting receiver is a WL extension of the minimum 
mean square error (MMSE) receiver and jointly mitigates multi-user interference and IQI by processing 
the real and the imaginary parts of the received signal separately. The IQI arising at both the base station 
(BS) and the user terminals (UTs) is then taken into account. The considered channel state information 
(CSI) acquisition model includes the effects of both estimation errors and pilot contamination, which 
is caused by the reuse of the same training sequences in neighboring cells. We apply results from 
random matrix theory to derive analytical expressions for the achievable sum rates of the proposed IQI 
aware and conventional IQI unaware receivers. Our simulation results show that the performance of the 
proposed IQI aware WLMMSE receiver in a system with IQI is close to that of the MMSE receiver in 
an ideal system without IQI. Moreover, our results for the sum rate of the IQI unaware MMSE receiver 
reveal that the performance loss due to IQI can be large and, if left unattended, does not vanish for 
large numbers of BS antennas. 
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I. Introduction 


M ULTIPLE-input multiple-output (MIMO) techniques have become a central part of 
modern communication systems such as Long Term Evolution (LTE) and WiMAX. With 
MIMO technology, high throughput and transmission reliability can be achieved. An emerging 
research field in wireless communications are so-called massive MIMO systems [HI, [l2l. Massive 
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MIMO systems employ a large number of antennas, e.g. one hundred antennas or more at the 
base station (BS), and ean aehieve very high speetral and energy effieieneies [l3]|. Moreover, 
in massive MIMO systems, the transmit power of the BS and the user terminals (UTs) ean be 
deereased by inereasing the number of antennas at the BS [|3||. These and other desirable features 
render massive MIMO a promising teehnology for future wireless eommunieation systems. In 
this paper, we eonsider the uplink of a multi-eell massive MIMO system. In an uplink single- 
eell massive MIMO system, whieh embodies a multiple aeeess ehannel (MAC), sueeessive 
interferenee eaneellation (SIC) is the optimal deteetion seheme [|4l. Nevertheless, sinee SIC is not 
feasible in most practical systems due to its high computational complexity, linear detectors such 
as matched filter (MF) and minimum mean square error (MMSE) detectors are often preferred 
as they provide a good trade-off between performance and complexity Ifll. If3l. 

However, hardware (HAV) impairments, which exist in all practical systems, can severely 
degrade the performance of linear detectors. One of the most important HAV impairments in 
digital communication systems is in-phase/quadrature-phase (FQ) imbalance, which arises in 
direct conversion transceivers [l5]|. In such systems, the real and imaginary parts of the received 
RF signal are mixed with the high-frequency carrier signal and its 90° phase-shifted version, 
respectively, to produce the baseband signal. Ideally, the phase difference between the carrier 
signal and its phase-shifted version is exactly 90° and the mixers for the real and the imaginary 
part have the same amplitude gain. However, in practical systems, both amplitude and phase 
mismatches between the real and imaginary parts occur, which leads to an I/Q imbalance (IQI) 
in each antenna branch at the base station and in the RF chain of each UT, and impair the 
received data vectors lIH, |171. Thus, for reliable detection, besides the mitigation of multi-user 
interference, IQI compensation is necessary as well. One approach to overcome the negative 
effects of IQI is to measure and compensate the individual IQIs in each antenna branch. However, 
this solution becomes very costly in massive MIMO systems, where the BS may be equipped 
with several hundred RF chains. Another approach, which we consider in this paper, is joint 
data detection and IQI mitigation. In this case, the equivalent channel, which comprises the 
actual channel and the IQI, is estimated at the BS. For this estimation, we exploit the received 
training sequence and the channel statistics which also include the effect of IQI. The detection 
matrix is then constructed based on this equivalent channel estimate and can be used for joint data 
detection and IQI mitigation. To the best of the authors’ knowledge, the problem of joint channel 
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estimation, data detection, and IQI mitigation for uplink multi-cell massive MIMO systems has 
not been investigated yet. 

Recently, a WLMMSE beamformer for systems with IQI has been proposed in |[8l. Here, 
the authors propose a beamforming scheme, where a multi-antenna receiver suffering from IQI 
detects signals coming from a specific direction while suppressing signals arriving from other 
directions. Another related work is lIH, where the authors propose a detection scheme for uplink 
single-cell multi-user-MIMO (MU-MIMO) systems impaired by IQI. In ifTOll . the authors model 
the residual HAV impairments, which remain after compensation, in massive MIMO systems as 
an additive Gaussian impairment and derive a capacity bound. Moreover, in our recent work 
m, we have proposed an IQI aware precoder for downlink single-cell massive MIMO systems 
assuming perfect CSI and IQI present only at the BS. Another recent work is |[T^ . where the 
authors consider the uplink of a single-cell massive MIMO system with IQI present only at the 
BS. In this paper, we propose a widely-linear MMSE (WEMMSE) receiver for CSI acquisition 
and data detection in an uplink multi-cell massive MIMO system, where both the BS and the 
UTs are impaired by IQI and both the received data and training signals for channel estimation 
are affected by IQI. WE filtering was introduced in ifT^ and is used to estimate complex signals 
by filtering the real and imaginary parts separately. Several works have employed WE filtering in 
single-user MIMO systems, cf. [[T4ll . [flSll . WE processing results in a higher performance than 
strictly linear processing, if rotationally variant signals are involved, which is the case when IQI 
is present. This motivates the use of WE filtering for channel estimation and data detection at 
the BS of uplink multi-cell massive MIMO systems suffering from IQI. 

In contrast to flUl, dH, and [fTTll . where a single-cell system with perfect CSI was considered, 
in this paper, a more sophisticated system model is adopted, which includes the effects of multi¬ 
cell interference and CSI imperfection originating from pilot contamination, channel estimation 
errors, and IQI. Eurthermore, contrary to ifT^ . where a single-cell system with ideal UTs was 
considered, our system model includes multi-cell interference and takes IQI at both the UTs and 
the BSs into account. Moreover, as opposed to ffTOll . where the residual HAV impairments were 
modelled by an additional additive Gaussian noise term and a non-augmented system model 
was employed, we adopt an augmented system model, which is essential for the analysis and 
mitigation of IQI. In addition, we use results from random matrix theory and provide analytical 
results for the sum rate performance of the proposed IQI aware (IQA) and conventional IQI 
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unaware (IQU) receivers. 

This paper is organized as follows. In Section |II1 the system model is presented. In Section 
Hm the conventional IQU-MMSE receiver is investigated and an analytical expression for its sum 
rate in the presence of IQI is derived. In Section |IVl we introduce the proposed IQA-WLMMSE 
receiver, and present an analytical expression for the corresponding sum rate. Numerical results 
are provided in Section |Vl and conclusions are drawn in Section |Vll 

Notation: Boldface lower and upper case letters represent column vectors and matrices, respec¬ 
tively. blockdiag (Q^^,..., Q^y) is a block diagonal matrix with matrices Q^,..., Q^y on its main 
diagonal. Ik denotes the K x K identity matrix and [A]^., [Aj.^, and [A]^; stand for the kth 
row, the /th column, and the element in the kth row and the /th column of matrix A, respectively. 
(•)* denotes the complex conjugate and det(-), H-IH, tr(-), (-j^, and (•)'^ are the determinant, 
Erobenius norm, trace, transpose, and Hermitian transpose of a matrix, respectively. and 
Q={-} denote the real and imaginary parts of a complex number, respectively. IE{-} stands for the 
expectation operator and CJ\f (u, ^>) denotes a circular symmetric complex Gaussian distribution 
with mean vector u and covariance matrix ^>. Moreover, a.s. stands for almost sure convergence. 

II. System Model 

In this paper, we consider the uplink of a multi-cell massive MIMO system with universal 
frequency reuse. The number of cells is denoted by L, and in each cell, K single-antenna UTs 
simultaneously transmit data to a BS with N antennas. K and N are assumed to be very large 
with their ratio /3 = K/N being constant. Eurthermore, we assume a block fading channel. 
The channel matrix between the UTs in the /th cell and the BS in the Ah cell is denoted 
by Hi,/ = [\i^i.. .\i^k] e Here, e is the channel vector 

between the kth UT in the /th cell and the BS in the Ah cell, where ~ CAA(0,lAr) 
and f, represents the channel covariance matrix. Since the 

detection schemes considered in this paper, i.e., IQU-MMSE and IQA-WEMMSE detection, 
have fundamentally different structures, we adopt two different representations for the system 
model, namely a complex-valued and a real-valued representation, which are presented in the 
following subsections. 
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A. Complex-Valued Representation 

In this subsection, the complex representation of the system model, which is used for the 
IQU-MMSE detector, is introduced. The transmitted data symbols of the K UTs in the /th cell 
are stacked into a vector, which is denoted by d; = ..., di^xY ~ (0) ^k), where di^k is 

the data symbol transmitted by the k\h UT in the /th cell. The received signal at the Ah BS can 
be modeled as 


Ti = 




A^i 


v« = l / \ 1 = 1 / 


( 1 ) 


where puL denotes the uplink transmit signal-to-noise ratio (SNR), and = diag(.^^_i^i,..., 
U,i,k) and Hb,i = dicLg{^B,i,i, ■ ■ ■ Ab,i,k) with ^A,i,k = cos (0z,fc/2) + sin (0i,fc/2) and 
^B,i,k = ^i,k cos [9i^k/‘2) — jsin(0; fc/2) representing the IQI at the fcth UT in the /th cell. ^ 
and ii^k denote the phase and amplitude imbalances at the corresponding UT, respectively. 
The IQI at the Ah BS is modelled by diagonal matrices ^A,i = diag • • •, ^A,i,Ar) and 

^B,i = diag('0B,i,i,... where ipA,i,n = cos(6'i,„/2) + je^^nsin (6'i,„/2) and iljB,i,n = 

€i^n COS {Oi^n/ 2) — j sin (^j „^/2) with 9i^n and being the phase and amplitude imbalances at 
the nth antenna branch of the Ah BS, respectively. If IQI is absent, we have = 0, 9i^n = 
0,V/ e L}, Vn e iV} and = 0, 9i^k = 0, V/ G {1,..., L} , V/c G iT}. 

Zj ~ CAA(0,lAr) represents the complex additive white Gaussian noise (AWGN) at the Ah BS. 

B. Augmented Real-Valued Representation 

Since IQI affects the real and imaginary parts of a signal differently, IQA channel estima¬ 
tion and data detection should allow for processing the real and imaginary parts of received 
signals individually. Hence, we use an augmented representation for the IQA system model, 
where the real and imaginary paAs of the signals are stacked together. More precisely, the real 
and imaginary parts of the independent and identically distributed (i.i.d.) zero-mean complex 
Gaussian transmit data symbols of the K UTs in the /th cell are stacked into the augmented 
vector dz = di^, ■ ■ ■ ,di, 2 K = dj]"^ G where d^^ = ■ ■ ■ Ari^kV and 

di; = [di^ , dl^ contain the real and imaginary paAs of the K complex transmit data 
symbols in the /th cell, respectively, and E |dzd; I — 0.5 l2K- The augmented data vector 
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received at the ith BS, fj G 


\,2K 


, from all UTs in the L cells can be expressed as 

= K + Zi , 


( 2 ) 


, i=i 


where r^. and ri. are the real and the imaginary parts of the received signal r, in ©, respectively. 


^2Nx2K 


with 


Here, the augmented real-valued channel matrix Hj; is given by 

Hi,, Hi 

Hr. ^ and Hr ^ being the real and imaginary parts of the complex-valued channel matrix Hr; = 


Hr 




Hi 


H 


R-i, 


j'Hi. respectively. The augmented real-valued vector Zj = [zr.Zr] G 


r,2N 


contains 

T 


the real and imaginary parts of the complex-valued AWGN vector z* = [ 2 :^ 1 ,..., ^rjv] = 
zr. + jzi. ~ CAA(0,lAr) at the i\h BS, respectively. G M 2 Afx 2 Ar models the 

IQI at the Ah BS. Here, is a permuted version of = blockdiag (Aj 1,..., Aj iv), where 
Arji, n G {1,..., A^}, represents the IQI of the nth RF branch of the Ah BS and is given by 


A • — 

■^i.n — 


The elements of permutation matrix 11 G M 2 iVx 2 Af defined as 

1 if p = 2n — 1, g = n, Vn G {1,..., A} , 

1 if p = 2n, g = n + A^, Vn G {1,..., A} , 
0 otherwise. 


(3) 


PI.,, = < 


(4) 


This permutation is required, since the stacked received data vector fj contains the in-phase 
and quadrature-phase components in its upper and lower parts, respectively. Moreover, in dH), 
H; = HHiH ^ G k 2 Rx 2 r (jenotes the stacked IQIs of the K UTs in the /th cell, where the 
permutation matrix H is similarly defined as H and is obtained by replacing N with K and n 
with k in (@1). Furthermore, H; = blockdiag (3;,i,..., Srr), where H; ^ G is the IQI matrix 
of the fcth UT in the Ah cell, which can be expressed as 

3 ? {^A,l,k + {^B,l,k — ^A,l,k} 

^ {^A,l,k + ^B,l,k} ^ {^A,l,k — ^B,l,k} 




( 5 ) 
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III. IQI Unaware MMSE Receiver 

In this section, as a performance benchmark, the sum rate of a conventional IQU-MMSE 
receiver comprising an IQU-MMSE channel estimator and an IQU-MMSE data detector is 
investigated in the presence of IQI. Eor IQU estimation and detection, we adopt the conventional 
MMSE estimator and detector, respectively, which are not designed for IQI mitigation. 

A. Channel Estimation 

In this subsection, channel estimation for an IQU system in the presence of IQI is presented. 
Eor channel estimation, at the beginning of every coherence interval, training sequences are 
transmitted by all UTs to their serving BS. Due to the limited length of the coherence interval, 
there are not enough orthogonal training sequences for all UTs in all cells. Hence, UTs with 
the same index in different cells use the same training sequence IfT^ . This leads to a corrupted 
channel estimate and this effect is known as pilot contamination in the massive MIMO literature 
m- Since we consider full pilot reuse, when pilot contamination is present, UTs having the same 
index k in different cells employ the same training sequence G where T is the length 

of the training sequence. The received training signal at each BS is multiplied by the original 
transmitted training sequence to eliminate the interference caused by other UTs. Considering 
O, and the orthonormality of the training sequences, i.e., = 1, x^Xj = 0, A; ^ j, we have 

the following expression for the received training sequence of the fcth UT in the Ah cell 

L L 

yi,k =YiXfc = ^A,i f y/PTRiii,l,k{^A,l,k + ^B,l,k) + Wj) + ^ (G,«,fc+ 

l=l 1=1 

where Yj G is the received training signal at the Ah BS. Here, ptr is the transmit training 

signal-to-noise ratio (SNR) and Wj = W^x^ ~ CAf (0, Itv), where Wj G is the AWGN at 
the Ah BS during the training period. In this paper, for IQU channel estimation, we assume that 
the estimator tries to estimate gi i^k = ^A,i,k'^A,ii^i,i,k as the desired channel between the kth UT 
in the Ah cell and the Ah BS. Since the IQU estimator does not process the real and imaginary 
parts of the received training sequence separately, it can consider only one component, i.e., 
^A,i,k^A,ii^i,i,k of the equivalent channel vector. We note that with the conventional complex¬ 
valued system model, which is assumed for the IQU-MMSE estimator, it is not possible to fully 
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model the equivalent channel vector, which comprises both the actual channel and IQI. Taking 
this into account and considering ^ as the observation, the MMSE channel estimate can be 
expressed as IfTTl 




= $ 


^i,i,ky i,k 


$ 


y i,ky i,k 


"y i,ki 


( 7 ) 


where $ 0 - , v , is the cross-correlation matrix of the desired channel estimate and the observa- 
tion, and given by 


$ 


Si.i.ky i,k 




A,i‘ 


( 8 ) 


The auto-correlation matrix of the received signal in ® can be expressed as 

L 

^yi,kyi,k = ® = X^PTR(|^A,Z,fc + + \^A,l,k + ^B,l,k\‘^^ 

1=1 

+ ^B,i^B,i. ( 9 ) 


X 


Now, we substitute ([8]) and dH) into dU) and obtain the following expression for IQU-MMSE 


estimation of the fcth UT’s channel vector 

L 


Si,i,k —^i,k[ y^hz,Z,fc('CA,Z,fc + ^B,l,k) + -f '^K,l,k{^A,l,k + CB,l,k) + 


1=1 


1=1 


:W,- 


\/PTR 

where deterministic matrix Oj ^ is given by 

^i,k =^A,i,k {^A,i,k + CB,i,k) ( (l^A,Z,fe + ^B,l,k\‘^'^A,iRi,l,k^A 


( 10 ) 


A,i 


1=1 


-1 


\Ul,k + ^B,l,k\^^B,i'Rll,k^B,^) + - 

Ptr 


( 11 ) 


If pilot contamination is absent, (fT^ reduces to 


Si,i,k =^i,k \ ^A,i[\i,k{^A,i,k+^B,i,k)- 


i/PTR 


-w, 


+ ^B,i ( {^A,i,k + ^B,i,k) ■ 


\/Ptr 


-w,- 


( 12 ) 


where the deterministic matrix ^ is equal to ^ if we set L = 1 and / = i in (fTTI) . 
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B. Data Detection 

In this subsection, we investigate IQU-MMSE data detection. The IQU-MMSE detector adopted 
here is the conventional single-cell MMSE detector. The IQU-MMSE detection vector for the 
/cth UT at the ith BS is given by 


^ H \ 

j -f -— In 
' Pul 


(13) 


where the fcth column of the estimated channel matrix Gj j is ga^ and given in (fTOl) and 
(fT^ for the cases with and without pilot contamination, respectively. Thus, the detected signal 
corresponding to the kth UT in the iih. cell at the output of the IQU-MMSE detector of the iih. 
BS can be expressed as 

=Ui,fcri = ^ -f -f Zi) + ^B,i f\/PULH*; (H^ ;d;* -f 


1=1 


■^BJ 


dz) + z* 


(14) 


C. Asymptotic Sum Rate Analysis 

The performance metric considered in this paper is the ergodic sum rate, which is a commonly 
used metric for performance evaluation of wireless communication systems. Eor the dh cell, the 
ergodic sum rate is given by ^ 

flW = y]E{log,(l + SINRg")}, (15) 


k=l 


where the expectation is taken with respect to the channel realizations. SINR]*^^ is the signal- 
to-noise-plus-interference ratio (SINK) for the Pth UT in the dh cell at the dh BS and given by 


sinrJ^^ = 




IQU 


i.k 


jlQU ylQV ’ 


(16) 


where Slf, , and Zlf are the useful signal power, interference power, and noise power for 

the Pth UT in the dh cell, respectively. In this paper, using results from random matrix theory, 
first, an analytical expression for SINR^*^^ , the asymptotic value of SINR]*^^ for large numbers 
of antennas N, is derived. Then, using SINR^^^ , the asymptotic sum rate is calculated as 


R 


K 

,IQU° 


= }_^log2(l + SINR 


IQU° 

i,k 


k=l 


(17) 
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In the following Theorem, we provide an analytieal expression for the asymptotie SINK of the 
IQU-MMSE deteetor. 

Theorem 1: In an uplink multi-eell massive MIMO system employing an IQU-MMSE reeeiver 
at the BS in the presenee of IQI, the asymptotie SINK eorresponding to the kth UT in the ith 
eell for N ^ oo is given by 


r-IQU° 

CT]^-plQU° _ _U;_ 

^IQU° ^IQU° ’ 


(18) 






where the asymptotie 
the asymptotie noise power, 

^iQU ^ Pul 

iV^oo “ (1 + 

Pul 


TOTT° TOTT° 

useful signal power, S-^ , the asymptotie interferenee power, , and 


-, , are given by 


CIQU“ _ j: oIQU _ _ 


1 3-A,2,/C‘^2,2,/c 


(19) 


iT = li- iT = 


N^oo (1 + 5. 


2 2 ^ ^ 

iB,i,k^\^i]k + CA,i,k^\,i,k + ) 

q=l 1=1 
{q,l)jt{k,i) 


1=1 


^ (Pilg + PmI) + + \^B,l,kf) (\7,lc + 

1=1 

= hm Z^T = 

i.k XT , i.k 


N^oo 


(1 + Si^kY 


tr ( 0 M.^ (r£^+r:;r) 


( 20 ) 


( 21 ) 


Here, <5^,^, and are defined as 


Si,k = —tr 


N 

1 


•i,i,k - i^A,i,k + ^B,i,k) tr 




( 22 ) 

(23) 

(24) 


'^i,l,q 


tr fc) 


23fJ< 


(1 + kY 

(R^T 9 ^A,ir-fca,fc^A,i) 


iV2 (1 + 5, 


(25) 




and g is identieal to ^ after replaeing in (l25l) superseript (A) by {B) and Ri,;,^ by Rj^^g. 
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Furthermore, for the case of pilot contamination, = E obtained as 


L 



1=1 



(26) 


and for the case without pilot contamination, we have 



(27) 


Moreover, F* is given by T in Lemma [3] in Appendix A for a = 1/(A^pul), 


and B = Otv- Furthermore, F' ^, and F^^^ are given by T' in Lemma |4] in Appendix A 
after setting A^. = &i^i^k/N and replacing C with ©*_*_*,, and respectively. 


Proof: Please refer to Appendix B. 


D. Asymptotic Sum Rate Analysis for Single-Cell Case 

Due to the very general setting considered in Theorem [H the obtained analytical expression 

for the asymptotic SINK of the IQU-MMSE detector is quite involved. Nevertheless, using these 

analytical results for performance evaluation is still much more convenient than performing 

lengthy Monte-Carlo simulations. However, to get some insight for system design, in this 

subsection, we provide analytical results for the simplified single-cell case with i.i.d. channel 

vectors and perfect CSI. In particular, we investigate the impact of the IQI at the BS and at the 

UTs separately to determine whether the IQI at the BS or at UTs is more harmful. 

Corollary 1: In an uplink single-cell massive MIMO system with i.i.d. channel vectors, perfect 

CSI, and IQI only at the BS, the asymptotic SINK of the kih. UT for the IQU-MMSE detector 

for K, N ^ oo, K N, and e„, On 1, is given by 



SINR° 


TQU-BS “ 


where e„ and On are the amplitude and phase imbalances of the RE chain of the nth antenna at 


the BS. 

Proof: Please refer to Appendix C. 
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Remark 1: For identical IQI at all BS antenna branches, i.e., e„ = e, = 0, Vn, the asymptotic 


SINK of the IQU-MMSE receiver for K, N ^ oo, K N, and e„, <C 1, is given by 

NpVL 


SINK- 


TQU-BS 


KpuL (e2 + i02)+i' 


(29) 


Proof: Substituting e„ = e,0n = 9,'in into (l2^ and considering e, 6> <C 1 leads to (l29l) . ■ 
Remark [H reveals that for a fixed number of users K, the SINK increases with increasing number 
of BS antennas N. 

Remark 2: In the absence of IQI, the asymptotic SINK of the IQU-MMSE receiver for 
K, N ^ oo, K N, is given by SINR^q-iqi = which is obtained by setting e = 6 = 0 

in (HH). 


Corollary 2: The asymptotic SINK loss of an IQU-MMSE receiver due to IQI in the uplink 
massive MIMO system defined in Corollary [H and for = e l,9n = 0 l,Vn e 
{1,.. .,N}, is given by 


lini A^inr = lim 


SINR^ 


No-IQI 


SINR°qu 


— -^PUL ( -f — 9 ^ ) -f 1. 


(30) 


Proof: Considering (|2^ . Remark [2l and performing simple mathematical manipulations 
yields (l30l) . ■ 

Erom (l30l) . it can be seen that the SINR loss of the IQU-MMSE detector compared to the ideal 
case without IQI does not vanish even in the asymptotic scenario where the number of the BS 
antennas is much larger than the number of the UTs. This motivates the need for a receiver, 
which mitigates both multi-user interference and IQI, cf. Section |IVl 

Corollary 3: In an uplink single-cell massive MIMO system with i.i.d. channel vectors, perfect 
CSI, IQI only at the UTs, and equal amplitude and phase mismatches, i.e., = e, 9k = 9,ik, 

the asymptotic SINR of the fcth UT for the IQU-MMSE detector for K, N ^ oo, K N, and 
e, 0 -C 1, is given by 


SINR' 


Np\]L 


IQU-UT 


iVpuL {P + \9^) + l 


(31) 


Proof: Please refer to appendix D. ■ 

Erom Corollary [3l we observe that, in an uplink massive MIMO system with IQI only at the UTs, 
for iV —)■ oo, the SINR does not increase with increasing number of BS antennas. Moreover, 
by comparing (|3TI) and (l29l) we conclude that assuming identical values for the amplitude and 
phase imbalances at the BS and the UTs, and K N, the SINR of a massive MIMO system 
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with IQI only at the BS is higher than that of a massive MIMO system with IQI only at the 
UTs. Thus, if not eompensated, IQI at the UTs is more harmful than IQI at the BS. 

IV. I/Q Imbalance Aware WLMMSE Receiver 

In this section, the proposed IQA-WLMMSE receiver, which comprises an IQA-WLMMSE 
channel estimator and an IQA-WLMMSE data detector, is presented. The IQA-WLMMSE 
channel estimator and detector process the real and imaginary parts of the received signals 
separately. The widely-linear filtering is necessary for mitigation of the IQI, since the real and 
imaginary parts of the signals are affected differently by the IQI. 

A. Channel Estimation 


The proposed IQA-WLMMSE channel estimation scheme is the widely-linear extension of 
the strictly-linear IQU-MMSE channel estimator introduced in Section IIII-AI and performs joint 
IQI compensation and MMSE channel estimation. Since we assume full pilot reuse, UTs in 
different cells with the same UT index k use the same training sequence resulting in the 

same augmented training sequence = 

system model. In order to mitigate the interference from other UTs, the received training signal, 
Yj, is multiplied by the augmented training sequence X^. Considering (jH), and taking into 
account the orthonormality of the pilot sequences, the training signal of the fcth UT received at 
the Ah BS is given by 

_ L 

^i,k = YjXfc = s/k^l,k + (32) 

l=l 


3?{xfc} -^{Xfc} 

^3{xfc} 3?{xfc} 


n2Tx2 


in the augmented 


with Wj = WjXfc, where Wj e £^‘iNx 2 T augmented AWGN matrix at the Ah BS. In (l32l) . 
Hj / fc £ represents the augmented channel between the kth UT in the /th cell and the Ah 

BS and is defined as 



3 ^ 












where ~ C-AA(0,lAr) and is the square root of the augmented channel covariance 
matrix between the kth UT in the Ith cell and the Ah BS. Here, G M2iVx2A given by 


R 


i,l,k — 




(34) 
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The proposed WLMMSE ehannel estimator estimates the equivalent augmented channel matrix 
of the kih. UT, which is given by Zi,i,k = Zi,i,k+K\ = ^i^i,i,k^i,k, where and 

Hj represent the IQI at the i\h BS, the actual augmented channel of the kih. UT, and the IQI at 
the kih UT, respectively. Note that for the IQ A-WLMMSE detector, both Zi,i,k Zi,i,k+K 
required and are used for detection of the real and imaginary parts of the data signal, respectively, 
c.f. Section IIV-BI The proposed WLMMSE channel estimate is given by 


§i,i,fc+ir 


= $- . . 


-1 


Y,- 


2,A:? 


(35) 


where $ 




is the auto-correlation matrix of the received training sequence Yj ^ and is given 


by 






i=i 

T „ _ ~ ~ T 


N—^oo 


Ptk 




(36) 


1=1 


Here, = E is obtained as 

"^1,1,k = ^i,l,k^ \'^i,l,k^l,k'^J,k^i,l,k\ ^J,l,k ./ ^ 0.51| || 


(37) 


where = 




mathematical operations. In (1351) . 


Gri i Yj t 


and the last step involved Lemma [Hand straightforward 
is the cross-correlation between the received training 


sequence and the desired channel estimate, and is given by 


Substituting (1381) and (1361) into 


estimate of the channel between the kih UT in the Zth cell and the Ah BS 

^i,i,k I 

\l=l 


Si,i,k Si,i,k+K 


(38) 

leads to the following expression for the IQA-WLMMSE 
e Zth cell and the Ah BS 

PTR^iili,l,k'^l,k + ^iWj I , 


(39) 
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where the IQA-WLMMSE channel estimator j ^ $ 




$ 


i,i,k ^ i,k \ '^i,k^i,k 


-1 


as 






yl'=l 


PTR 


can be expressed 


(40) 


Remark 3: One of the main features of the proposed IQA-WLMMSE channel estimator is 
that it does not require the explicit knowledge of the IQI at the BS and UTs; indeed, only 
the equivalent channel covariance matrix and the autocorrelation matrix of the received training 
signal are needed. This can be observed from (l39l) . where the augmented channel estimate 
is obtained by filtering the received training signal, i.e., the term in parenthesis, with the 
IQA-WLMMSE estimator The IQA-WLMMSE estimator comprises two components: 

= and = = . The first component = , which is given in (l3^ . can be rewritten 

as PtrE {Si,i,kSjik}^ which is the product of the training SNR and the covariance matrix of the 
augmented equivalent channel and both can be estimated using SNR estimation techniques 
and channel statistics estimation techniques, respectively. The second component in the IQA- 
WLMMSE channel estimator is , which is the autocorrelation of the received training 

signal. 

Eor the case without pilot contamination, the IQA-WLMMSE channel estimates are given in 
(1^ after setting L = 1 and I = I' = i by (1^ and (l40l) . 


B. Data Detection 


The proposed IQA-WLMMSE data detector is the widely-linear extension of the conventional 
MMSE detector considered in Section IIII-R[ and employs the estimate of the equivalent channel. 


G,,, = 


• • • 1 


, which comprises the actual channel, the IQI at the BS, and the IQI 
at the UT. The IQA-WLMMSE detector includes the filter vectors for the real and imaginary 
parts of the signal of the fcth UT at the fth BS and is given by 


^i,k+K 


^i,i,k 

^,i,k+K 




(41) 


where gj ^ ^ and gi^i^k+K the kib and the {k + K)ib columns of the estimated augmented 
channel matrix respectively, and are given in (l39l) . Hence, the decision variable at the output 
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of the IQA-WLMMSE detector at the ith BS corresponding to the fcth UT can be expressed as 


JIQA 




\/puL^iHj 


(42) 


, ;=i 


. z=i 


C. Asymptotic Sum Rate Analysis 


The ergodic sum rate of the IQA-WLMMSE receiver is given by 

K 


fl|«* = y^E{log2(l + SINR:«'')}, 


(43) 


k=l 


where the expectation is taken with respect to channel realizations. Here, is the SINK 


of the fcth UT in the ith cell at the ith BS and is defined as 

oIQA 

sinr!*?^ = 




IQA ^IQA ’ 
i.k 




(44) 


where and are the useful signal power, interference power, and noise power 

of the kth UT at the Ah BS, respectively. Using again results from random matrix theory, we 
will show that the asymptotic sum rate of the IQA-WLMMSE detector can be expressed as 

^IQA° ^ ^ ^ SINR|Q/°) , (45) 


k=l 


where the asymptotic SINK expression SINR]^^ is provided in the following theorem. In Section 
rvl we show that the derived asymptotic sum rate accurately predicts the ergodic sum rate, which 
is obtained through lengthy Monte-Carlo simulations. 

Theorem 2: In an uplink multi-cell massive MEMO system employing an IQA-WLMMSE 


receiver, for N ^ oo, the asymptotic SINK of the fcth UT in the Ah cell is given by 


SINRg^° 


r-IQA° 

‘^i,k 


rIQA° I ^IQA° ’ 


(46) 


^i,k 


where the asymptotic useful signal power, , the asymptotic interference power, , and 


^i,k 

the asymptotic noise power, , are given by 

( 

\Xi,k\ , lu 

O, = U.DDTIT, — 

N^oo 


*^3^° = = 0-^^uL 




V 


1 + 


+ 


i,k+K\ 


1 -f (5, 


i,k+K 


(47) 
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ilT’ = J™,^2* = 




2 


^i,l,q 



\ l = I 

q^{k,k+K} 



, ^i,l,q'^i,i,q' 

1 + ^i,fc 


= lim = 0.5 


/ ~ ( 2 ) 

Here, we use Xi,k = tr ( ^ 




where and are defined as 

L 


(i + *.,) 

iV / 
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~T 


^ z 

1 f~ ^i,i,q' 

' 7 1 ‘"2^,9' + -7 „ N 2 

*=1 '?'=! ( 1 + '5i,k ) \ ( ^ 2 ” ) 

q>^k,2k ^ ^ 

2 

^i,l,q' ^i,i,q' 


1 f -^*-29' ^i,i,q' 

7 I ^ I 7 ~ ^ 

( 1 + Si,k+K) V f 1 + 5i,fc+A: 



tr(i.2,f'",) 


] 

^4iV2 

1 + Si^k^ 

2 ' 

4Ar2 

+ Si^k+K ^ 

2 

/ 


(49) 



’ ^i,i,q 
+ - 


^ tr 





2iVpuL 


and 


■l2Ar 


(50) 




H 1 ~ ~H 


^tlk - 


i,i,k ^0-5 ^ ^ - 

0.5 ^[ 1 ^ 2 ,/c]]^ ~t~ [‘^^,A:]2 1 ^ 


^l.,k ( 51 ) 


\l^) S72,i,fc^2Ri,i,fc^2^. (52) 

+ [^i,k]li with [Si,k]l2 + 

1521) . respeetively. In (l50l) . 


Furthermore, ^nd obtained by replaeing [3/,^]^ ^ + [^i,k]l^i with [1 

[3z,fc]2 2 and with + [^iA \2 lEB and ([521), respeetively. 

Si k and Si }i-\-K tVn=» colntionc to tVio follo\x/infT flvoH-ooir*^ 


diiu ^i^i^k+K uuLdiiicu uy icpidcmg ^ [=/,/cJ2^i vviiu 

+ [32,fc]2^^ with [^ix ]\^2 + [“i,fc] 2,2 lEB and ([521), respeetive 
the solutions to the following fixed-point equations 



(53) 


(54 


given by tr ^Bf/(4iV^), where B is equal 
I f ^ g is given by T' in Lemma |4] after replaeing 
Ai. = $!V,./(2iV). Similarlv. C ',.. RZ and 6', 


Moreover, ^i^i^q, Ri,i,q, and ^i^i^q in (l48l) are given oy tr /(^4iv‘j, wnere ±3 is equs 

Rj j g, and fli^k^i,i,q, respeetively, and f- ^ is given by T' in Lemma |4] after replaein; 

by 2N and setting T = f,, C = and = ^fJ^k/i‘^N). Similarly, and ^ 

are given by tr /(4iV^), where B is equal to Ri,i,g, ^^i^lk^ and fli^k^i,i,q, respeetb 
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and is equal to T' in Lemma 0] after replaeing N by 2N and setting T = Fj, C = 
and Afc = ^^Hk+K/ (2-^)- Furthermore, T'”g in (1^ is equal to T' in Lemma |4] after replaeing N 
by 2N and setting T = f C = and (2-^)- We note that if pilot eontamination 

is not present, I and L in (ISTI) are set to i and 1, respeetively, and is obtained from (l40l) 
after setting I' = i and L = 1. 

Proof: Please refer to Appendix E. ■ 


D. Asymptotic Sum Rate Analysis for the Single-Cell Case 


Although the provided asymptotie sum rate expression is easy to evaluate numerieally, sinee 

Theorem [2] eonsiders a very general ease, it does not offer mueh insight for system design. 

Henee, in order to get some insight regarding the influenee of IQI on the performanee of uplink 

massive MIMO systems employing the IQA-WLMMSE reeeiver, similar to the analysis for the 

eonventional IQU-MMSE reeeiver, we eonsider the simpler single-eell ease with perfeet CSI, 

and i.i.d. ehannels for all UTs. In partieular, in order to investigate the influenee of the IQI on the 

performanee, the eases with IQI present only at the BS and only at UTs are analyzed separately 

and their asymptotie SINRs and the eorresponding improvements eompared to the eonventional 

IQU-MMSE reeeivers are evaluated in the following Corollaries. 

Corollary 4: In an uplink single-eell massive MIMO system with i.i.d. ehannel veetors, perfeet 

CSI, and IQI only at the BS, the asymptotie SINK of the fcth UT for the proposed IQA-WEMMSE 

reeeiver for iC, —)> oo, K N, and e„, 1, is given by 


SINK 


IQA-BS° 

k 


a'pul (i + ^ Eti 

1 + jf En=i ((6 - 2 sin^ el + sin^ 9^) ’ 


(55) 


where e„ and 6n are the amplitude and phase imbalanees of the RE ehain of the nth antenna at 
the BS. 

Proof: Please refer to Appendix E ■ 

Remark 4: Eor the system deseribed in Corollary |4] and identieal IQI at all BS antenna 

branehes, i.e., en = e,6n = 0,\/n, the asymptotie SINR is given by 

IQA-BS° _ (1 + 


SINK 


(56) 


1 -f (6 — 2 sin^ P) -f sin^ 6 

Substituting e = 9 = 0 into (l5^ leads to SINR^*^^ = A^puL for an ideal system without IQI. 
Erom (15^ . we observe that with inereasing number of BS antennas, the sum rate of the proposed 
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IQA-WLMMSE receiver increases too. In particular, it can be shown that for e„, 1, which 

is valid for typical IQI values, the asymptotic SINK in (l5^ is smaller than but very close to 
TVpuLj i-e -7 the sum rate of the ideal system without IQI. In addition, comparing the SINK of the 
IQA-WLMMSE receiver with IQI present only at the BS given in (15^ and the corresponding 
SINK of the conventional IQU-MMSE receiver given in (l29l) . the following asymptotic SINK 
loss can be obtained 


lim A 

/3^0 


SINR-BS 


= lim 


_ (1 + eY (KpuL (e' + + l) 


(57) 


-^^0 SINR^^^ 1 + (6 — 2 sin^ 6*) + sin^ 0 

Erom (1571) . it can be observed that for systems with IQI only at the BS, the SINK loss increases 

with increasing SNR and increasing number of UTs. 

Corollary 5: In an uplink single-cell massive MIMO system with i.i.d. channel vectors, perfect 

CSI, IQI only at the UTs, and equal amplitude and phase mismatches, i.e., = e,9k = 0,V/c, 

the asymptotic SINK of the /cth UT for the proposed IQA-WLMMSE receiver for K, N 

oo, K N, and e, 0 <C 1, is given by 


SINR 


iqa-ut^ -^Pul (l + 2e^ (l 2 cos^ ff ^ 


2 (1 + e2) 


^PuL (1 - 2e) 

1 + iVpuL (1 + 2e)J ' V1 + iVpuL (1 - 2e) 


(58) 


Proof: Please refer to appendix G. ■ 

Substituting typical values for e and 6 into (15^ . it can be observed that similar to the system with 
IQI only at the BS, the asymptotic sum rate of the system with IQI only at the UTs increases 
with increasing number of BS antennas, and is smaller than but almost identical to the sum 
rate of an ideal system without IQI for e, P 1 and puL ^ 1- Moreover, considering (15^ and 
(|3TI) . we obtain the following asymptotic SINR loss of the conventional IQU-MMSE receiver 
compared to the IQA-WLMMSE receiver in a system, where the IQI is only present at the UTs 


SINR^Q^-^^° iVpuL (6^ + \e^) (( iVpuL (1 + 21) V / iVpuL (1 - 26) \ 
SINRfc^^“^^° 2(l + e2) l^Vl + ^PuL(l + 2e)y^Vl +ApuL(l-2e)y j' 

(59) 

Erom (l59l) . we observe that, if IQI is present only at the UTs, the SINR loss of the conventional 
IQU-MMSE receiver compared to the IQA-WLMMSE receiver increases with increasing number 
of BS antennas, N. Substituting typical values for e and 6, it can also be observed that the SINR 
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loss of the conventional IQU-MMSE receiver compared to IQA-WLMMSE receiver in systems 
with IQI only at the UTs is slightly larger than the corresponding loss in systems, where the 
IQI is present only at the BS. We validate this observation in Section |V] for multi-cell systems 
and more general settings. 


V. Numerical Results 

In order to evaluate the performance of the proposed IQA-WEMMSE receiver and to validate 
our analytical results, Monte-Carlo simulations have been performed. Here, we assume a system 
consisting of seven hexagonal cells with a normalized cell radius of one. Without loss of 
generality, we further assume that the central cell is the target cell. In each cell, there is a BS in 
the cell center and there are K UTs, which are uniformly distributed on a circle with a radius 
of 2/3. The channel model used here comprises path-loss, antenna correlation, and Rayleigh 
fading. Moreover, we assume that the BS employs a uniform linear array (UEA) and adopt the 
UEA channel correlation model used in ifT^ . In particular, we have 0Afx(iV-M)]> 

where Ci^i^k is the distance between the /cth UT in the /th cell and the Ah BS, and 0nx{n-m) and 
M are an N X {N — M) all-zero matrix and the number of dimensions of the antenna’s physical 
model, respectively. We adopt B = [b (0i),..., b (0Ar)], where the steering vector b (0^) is 
defined as b (0m) = [l, • • •, /\/M with 0m = —ti/2 + (m — 1) tt/M being 

the mth angle of arrival (AoA), and 7 and A being the antenna spacing and the wavelength, 
respectively ifT^ . 

In Eig. [H the ergodic sum rates of the IQA-WEMMSE receiver, the IQU-MMSE receiver, 
the MMSE receiver in the absence of IQI, and the MMSE receiver in an ideal system with 
perfect CSI are depicted. Here, we assume puL = 15 dB and pxn = 10 dB, and we further 
assume that the IQI is present at both the BS and the UTs. Moreover, except for the perfect 
CSI case, full pilot contamination is assumed. The number of UTs is set to iC = 10 and the 
amplitude and phase mismatches at the UTs and the different antenna branches of the BSs are 
randomly and uniformly distributed in the range of 0.15 < ei,ki^i,n < 0.2, 1° < Oi^kiQi^n < 
2°, V/c e {1, • • • , /T} , / e {1, • • • , , / e {1, • • • , L}, respectively. As can be observed from 

Eig. [H even small amplitude and phase mismatches lead to a high sum rate loss of the IQU- 
MMSE receiver compared to the ideal system without IQI. As expected from the analysis of the 
simplified single-cell channel model in Section IIII-Dl the rate loss associated with IQI does not 
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Fig. 1. Sum rate vs. number of BS antennas N for K = 10, puL = 15 dB, pxR = 10 dB, and pilot contamination 
with IQI present at the BSs and UTs. 

vanish even if the number of BS antennas is much larger than the number of UTs. For example, 
for N = 80, the rate loss of the IQU-MMSE receiver compared to the system without IQI 
is approximately 60%. Furthermore, as expected from the analysis of the simplified single-cell 
channel model in Section IIV-D[ the proposed IQA-WLMMSE receiver achieves a substantially 
higher sum rate than the IQU-MMSE receiver and closely approaches the sum rate of the MMSE 
receiver in an ideal system without IQI. In Eig. [H we also present analytical results for the 
asymptotic sum rates of the IQU-MMSE receiver and the IQA-WEMMSE receiver given in (fTTI) 
and (l45l) . respectively. Eor large N, a perfect match between analytical and simulation results 
is observed for all receivers. Nevertheless, even for small numbers of BS antennas, the match 
between asymptotic and simulation results is good. 

In Eig. [2l we investigate whether IQI at the UTs or IQI at the BSs is more harmful. To do so, 
we compare the sum rate performance of systems with IQI only at the BS (BS-IQI), IQI only 
at the UTs (UT-IQI), and IQI at both BSs and UTs (BSUT-IQI). Eor clarity of presentations, 
only analytical results are shown in Eig. [2l However, all results were verified by simulations. 
Here, we consider a system without pilot contamination but with channel estimation errors, and 
Pul = 15 dB and ptr = 10 dB. The amplitude and phase mismatches are generated in the same 
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Fig. 2. Sum rate vs. number of BS antennas N for K = 10, puL = 15 dB, ptr = 10 dB, and no pilot contamination. 

manner as for Fig. [T] From Fig. [2l we observe that if the IQU-MMSE reeeiver is employed, 
the system with IQI at both UTs and BSs yields the lowest sum rate, as expeeted. Furthermore, 
the system with IQI only at the BSs aehieves a higher sum rate than the system with IQI only 
at the UTs. We note that this effeet eould also be observed in Seetion IIII-Dl where analytieal 
expressions for the asymptotie SINRs in the simplified single-eell system were derived. In faet, 
the sum rate of the system with IQI both at the BSs and the UTs approaehes the sum rate of 
the system with IQI only the UTs for large numbers of BS antennas. A similar behavior ean be 
observed for the sum rate performanee of the IQA-WLMMSE reeeiver. Again, the system with 
IQI only at the BSs aehieves the highest sum rate followed by the system with IQI only at the 
UTs and the system with IQI both at the UTs and the BSs. We note that this behavior supports 
the results in [fT^ . where the authors elaim that in the asymptotie regime, where the number of 
BS antennas is very large, HAV imperfeetions at the UTs are more harmful than those at the 
BS. 

In Eig. |3l the influenee of the amplitude mismateh e on the sum rate of the eonventional IQU- 
MMSE and the proposed IQA-WEMMSE reeeivers is investigated. Eor eonvenienee, we assume 
that all UTs and BS antenna branehes have the same phase and amplitude mismatehes, i.e., 

h,k = = e, Kk = di,n = 9 = 2°, V/c e {!,■■■ ,K} ,i E {!,■■■ ,N}, I E {!,■■■ ,L}. The 
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Fig. 3. Sum rate vs. e for N = 100, K = 10, puL = 15 dB, pxR = 10 dB, and pilot contamination. 



Fig. 4. Sum rate vs. 0 for N = 100, K = 10, puL = 15 dB, ptr = 10 dB, and pilot contamination, 
number of BS antennas and the number of UTs are set to iV = 100 and K = 10, respectively. 
Moreover, we assume full pilot contamination and we further assume that the transmit data and 
training SNRs are puL = 15 dB and ptr = 10 dB, respectively. From Fig. [3l it can be seen that 
the sum rate of the IQU-MMSE receiver rapidly decreases with increasing amplitude mismatch 
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e. The proposed IQA-WLMMSE reeeiver performs signifieantly better than the IQU-MMSE 
reeeiver and its performanee loss eompared to the ideal system without IQI is small even for 
large amplitude mismatehes. 

The impaet of phase mismateh on the sum rate performanee is depieted in Eig. HI The same 
simulation parameters as for Eig. [3] are adopted. However, now the amplitude mismateh at all 
UTs and BS antenna branehes is set to e = 0.02, and the performanee is evaluated for different 
values of phase mismateh 6. Similar to Eig. [3l it ean be observed that the sum rate of the 
IQU-MMSE reeeiver deereases rapidly with inereasing mismateh. Moreover, the proposed IQA- 
WEMMSE reeeiver has a substantially higher performanee than the IQU-MMSE reeeiver and 
its performanee loss eompared to the ideal system without IQI is negligible even for large phase 
mismatehes. 


VI. Conclusion 

We have proposed an IQA-WEMMSE reeeiver for uplink multi-eell massive MIMO systems 
suffering from IQI at both the BS and the UTs. In the eonsidered system, CSI aequisition and 
data deteetion are affeeted by IQI, pilot eontamination, and multi-eell interferenee. The proposed 
reeeiver eomprises a ehannel estimator and a data deteetor and proeesses the real and imaginary 
parts of the reeeived signal separately. Our simulation and analytieal results show that, if left 
unattended, IQI eauses severe performanee losses even if the number of BS antennas is mueh 
larger than the number of UTs. The proposed IQA-WEMMSE reeeiver yields a substantially 
higher sum rate than the IQU-MMSE reeeiver and approaehes the performanee of the MMSE 
reeeiver in an ideal system without IQI. Eurthermore, we observed that for the eonventional 
IQU-MMSE reeeiver, IQI at the UTs is more harmful than IQI at the BS. We validated our 
simulation results by providing analytieal expressions for the asymptotie sum rate using tools 
from random matrix theory. 


Appendix A - Some Useeul Eemmas 


Lemma 1 ( / I79l Theorem 7]): Eet p, q G have mutually independent, i.i.d. zero-mean 

unit varianee Gaussian distributed elements and A G be a Hermitian matrix with bounded 

speetral norm, whose elements are independent of p and q. Then, 



N—>-oo N 


-tr (A) 


(60) 
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J- Ha r\ 

—p Aq-;■ 0. 

N W—>00 


(61) 


Lemma 2 ( I^Ui Rank-1 Perturbation Lemma]): Let A G q G and B G 

be a Hermitian nonnegative definite matrix. Moreover, let a and a' be positive real numbers. 
Then, 


tr ^A (B + alAf) ^ — A (B + a'qq'^ + alAr) 


-1 


< 


a 


(62) 


Lemma 3 ( / |Z71 Theorem 1]): Let A G and B G be Hermitian nonnegative 

definite and G G have random column vectors ~ CM {0, A^). Moreover, A, B, and 

Afc, Vfc = 1,..., A, have bounded spectral norms. Then, for positive a. 


-tr (a (GG^^ + B + aljv) 


-1 


N^oo N 


» 4tr (AT), 


where T is given by 


K 


-1 




iv^i + 4 


(63) 


(64) 


with 6k being the solution to the fixed-point equation 

K 


6k = —tr I A/i 


N^l + 6j 

j=i 4 


B + oT 


■N 


(65) 


Lemma 4 ( /EH/)-' Let A G B G and G G be defined as in Lemma 3. 

Then, for Hermitian nonnegative definite C G with bounded spectral norm. 


N 


where T' is defined as 


tr (A (GG^^ + B + alTv) ^ C (GG^ + B + al^) ^ 


TV—>-oo N 


> —tr(Ar) 


K 


T' = TCT + V 


( 66 ) 


(67) 


where T and 6k are given by (1^ and (1^ . respectively and 6' = [(i(,..., 6'^^] = (Jk — Y) x, 
where the elements of Y G and x G are given by 

Ivtr (AfcTAgT) ^ ^tr(AfcTCT). 

L AT (1 +5^)2 L Jfc ^ V fe ; 


( 68 ) 
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Appendix B - Proof of Theorem 1 

According to (fT4l) . for the IQU-MMSE receiver, the useful signal of the fcth UT reeeived at 
the ith BS ean be expressed as 

di,k = (^ H- In) ^A,ih.i,i,k + \/P^^^*B,i,kS^,i,k - ^n) 

\ Oiu. Oiu 


-1 


X ^B,i^i^i^k ) di,k- 


(69) 


Thus, in the large system limit for N ^ oo, the useful signal power of the kth UT reeeived at 
the Ah BS converges to 


= lim Pul 


N^oo' 
1 


(—^ + ^27—liv) ^A,iKi,k 


NpuL 


In) 


N ^ N 
2 '' 


NpvL 


+ 




N 


N 


(70) 


Applying the matrix inversion lemma If^ . Lemmas [Il[2l and[3l the first term on the right hand 
side of dVOl) ean be rewritten as 

C . . A-H 1 2 iA,i,kii,i,k A “llTt t. 2 C . . . \ 

^A,i,kbi^i^k . J- t I.t, u N ^i,k ^A,i^i,i,k a.s. 

I--^777- ^ 


N 


A^Pul' 


1 + jfS\',i,kd^i,k 


- 1 ; 


N^oo 


1 + 5- 


'i,k 


where Aj ^ is defined as 


1 - - H 1 1 

Ai.fc = - J^k^,^,kki,i,k + 


6i^k is given in (I22l) . and is given in (12^ for I = i. Similarly, the second term on the right 
hand side of (TTOl) ean be expressed as 




(71) 


(72) 


1 .. / 1 - ..H 1 




N 


-1 


-)■ 

N^oo 


\(-S) 


1 + 5- 


i^k 


(73) 


where X\^\ is given in (1241) for I = i. Next, we derive an analytieal expression for the asymptotie 
value of the interferenee power for N ^ oo. According to (fT4l) . the interferenee part of the 
reeeived signal of the fcth UT at the Ah BS is given by 


L K 


1=1 q=l 
{l,q)jl(i,k) 


N 


, -/ ^-4,«ijgj j fc/Gj^jGj j 1 

- ' -T^ + ^7-1 


A 27puL 


-1 


^ A,i^i,l,q~^ 


^B,LqSi,i,k ( Gi,iG 




N Npui^' 


-I 


N 


-1 
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X 


L K 

1=1 q=l 


^B,l,qSi,i,k f 


N 


N iVpuL 


^ T ^ \I> h -I- '^^-*-9^*-*-^ 

^A,W,l,q + 


N 


X 


GmG-,^ 1 


-1 


-1^1 ^B,iKi,g]dlq- 


N A^Pul 

In the asymptotic regime, when N ^ oo, the interference power in (1741) converges to 
^2"" = 


(74) 


1=1 q=l 
{l,q)^{i,k) 


( 

^A,l,qgi^i^k /GijiGj^j 


2 

L 


N N 




C* A-H P . .P 


N 


N 


liV \-l 


A^Pul 




+ 


^B,l,qSi,i,k Iv 


N 




i,l,q 


+ 


Iat \-l 


A^Pul 




') 


+PUL 


N NpuL 

'bS,*,fcSi i fc / Gi,iG j j Ijy \—1-j- -i 


c * P ■ P*^ 


N 


N 


N 


N A^Pul ' 




N 


X 


( GyG,, In x-i 


N A^pul ' 

The first term on the right hand side of (1751) can be rewritten as llT^ 


(75) 


L K ^ -H A. 


EE| 

1=1 q=l 
{l,q)^{i,k) 


L K 

1=1 q=l 
q^k 


^ A , l , qSi , i,k , liV \-1 t, 1 

' ' ) ^A,i\l,g\=} , 


L ^ ^ ~H 


N 


^A,l,kSi,i,k , liV ^-l.n 


N A^pul ' 


1=1 

l^i 


N 


^ N ^ATpuL^ *’ 




hbAj^igSi i fc / Gj^jGj j u 

■(- ^ + w:—) ^A,ihi,l,q 


N 


N A^Pul ' 


(76) 


Applying the matrix inversion Lemma ll22ll and Lemmas HI lH and lH yields 

^A,l,kg^,i,kfGi^iGii 1 a.s. ^A,l,k{CA,l,k +CB,l,k)^T^ {^A,i'Ri,l,k^%i^i,kTi) 


N 
1 + di^k 


N ^ iVpuL' 


7.^77- 

N^oo 


A^ (1 + 6i^k) 


(77) 


where and are given in (l22l) and (l23l) . respectively. On the other hand, for the second 
term on the right hand side of (l76l) . we define = E {gi^i^kk^,i,k}’ ^PP^Y the matrix inversion 
Lemma ll22l and Lemmas HI l2l lH and lH and obtain iH^ 


j\^^A,l,qgi^i^T 


N 




A,irii,/,q 


(1+ jfsZ,k\isi,i,ky 
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\^A,l,q\ 

|2 

(^Z,q^%^i,k,q®i,i,k\k,q^A,ii^i,l,q 

r,H *-i - 

^H,l,q ^ A,i'^ ^i,k,qoi,i,q \ 

^i,i,q'^'i,k,q^^^-)^i^‘^'i,k,q^i,i,q 

(1 + ()i, 


1 A2 ^ 

(1 + kgf 


- 23fJ< 


{sZAilq^A,iKi,q) 


A 


(A) 

i,l,q 


K 


1,1,q 


- 2^< 


x(A) ,(A) 
^i,l,qTi,i,q 

1 + & 


(1 + <5,,) 

It |2JA) 
|SA,Z,(j| Uilq 


a.s. ^ I^A,Z,ij| 


N^oo (1 + Si^kY 


AA) 

^i,l,q' 


i,q 


(1 + ^i,k) 


2 ’ 


(1 + ^i,Y 

where is given in (1^ and h.i^k,q is defined as 

1 ^ ’ 1 


^i,k,q ^Gi,iGi^i j^^Si,i,kSi,i,k Ar^A,gSi,i,q^ 


N 


iV' 


NpuL 


and 6i^q and A-f ^ are obtained by replaeing k with q in (l22l) and (l2^ . respeetively. In 
Ki^i^q, and are defined as 


^ N^oo 


tr (Ru,,>Plir:.t>t'/.,i) = C'A 




K. 


1,1,q 


— Kl,g'^li\lq®A,k\lqgi,i,q 


^A,l,q 




(78) 


(79) 


AA) 


(80) 

(81) 

(82) 


^2 t,k,q-^,^,>^ ^,k,qr,^,^,q ^2 

where F' ^ is given by T' in Lemma |4] with A, B, T, C, and being set to Qi^k^A,iRi,i,q'^A,v 
0, Fj, and Gi^i^k/N, respeetively, where 0 is an all-zero N x N matrix. Considering 

(iT^ - dSH) . and performing similar mathematieal operations for the remaining terms in (1751) . we 
obtain the expression in (l20l) for the interferenee power of the fcth UT reeeived at the ith BS. 
Furthermore, applying the matrix inversion lemma lf22l . the asymptotie noise power at the output 
of the IQU-MMSE deteetor eorresponding to the fcth UT in the fih eell is obtained as 


r^lQV° 

^i,k 


= lim 

N^OO 


sY,i,k\l^A,i^A,i\kki,i,k , sZ,k\l^B,i^B,i\lhi, 


+ 


.2 • 


(83) 


( 1 + jfk^,i,k\iUi,k) ( 1 + jfsY,i,k\isi,i,ky 

Then applying Lemmas 1-4, the asymptotie expression for the noise power at the output of the 
IQU-MMSE deteetor in (|2T]) is obtained from ([83]). This eompletes the proof. 

Appendix C - Proof of Coroffary 1 


Applying the matrix inversion lemma [22], and Eemmas 1, 2, and 3, the useful signal power 
of the fcth UT for IQU-MMSE deteetor in a single-oell uplink massive MIMO system with i.i.d. 
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channels, IQI only at the BS and perfect CSI (i.e., with ^ ^ A^i,i,k) is obtained as 


= 1™ Pul 


N—^OQ 


N 




Pul 


.h,- 


Pul<5' 


N\i,k^ A^i,k^ A^i,i,k 

~ (1 + ^ (1 + i) 


2 ’ 


where A^ = Ai,^ = A^ - jjgi^i^kSi,i,k^ and 5 is defined as 

^ = i™ ^Ki,k^AK,k^AKi,k = lim ^tr , 

N^oo I\ ’ N^oo I\ ^ ' 


and r is given by 


r = 


I3^a^ 
1 + 5 


" 1 


+ 


ApuL 




-1 


(84) 


(85) 


Assuming /) -C 1, 5 ^ 1, and applying Taylor series expansion, (1851) ean be rewritten as 

iVpuL/3*.*^^ ,3,^ 


r ApuL ivf — 


1 + 5 


Substituting (l84l) into (l8^ and eonsidering /) <C 1 and 5 S> 1, the following quadratic equation 
is obtained 

5^ + (1 - ApuLp) 5 + (aVul/^p' - Apuup) = 0, (87) 


where p = tr (^'a^a) /A and p' = tr ^(^'a^a)^^ /A. Solving (1871) and assuming /) —>■ 0 
leads to 

5 = (p + Vp2 + d/^p') -^ ApuuP- (88) 

Now, substituting (1^ into (l8^ and eonsidering (3 <^1 yields 


/3^-0 


■> ApulIat- 


(89) 


Furthermore, the interferenee power of the Pth UT in the eonsidered single-cell uplink massive 


MIMO system ean be expressed as 

^ 1 
rIQU“ _ 'Sp ^ ^ 

^i,k - 2^ Pul 

q=l^q^k 


K 


N 




5Z^ul 

9=1 


4h" (90) 


Using the matrix inversion lemma ll22l . and Lemmas [Il[2l and[3l the first term on the right hand 


side of (1901) ean be expressed as 
1 

iV ’ ’ 


K 


Pul 

g=l,g^fc 




AK^,qKi,q^^A^i 


q=l,q^k 


a.s. ^ /5 PUl5' 

(1 + 5)4' 


( 91 ) 
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Here, 5' is given by 




'{i+sf n'- 


Ltr 


/3-i>0 


> /x'iVV^jL, (92) 


where we used Lemma 4 and (I891) . Following a similar proeedure, and eonsidering <5 3> 1, the 


seeond term on the right hand side of (1901) ean be approximated as 

K ^ ^ K ^ 


q=l 


— h,^,- 4 A,- ^ bK i n ' 

TV ’ ’ ’ 




9=1 


a.s. ^ 

(1 + 5)2’ 

where 5" is given by 
lim — 

N^OO N 


(93) 




(1+5)’'A' 




/3^0 


(94) 


with +" = tr /N and we used Lemma 4 and (l89l) . Moreover, eonsidering (fT4l) 

for L = 1, perfeet CSI, and no IQI at the UTs, applying the matrix inversion lemma ^22^ and 


Lemmas [H [2l and |3l the reeeived noise power eorresponding to the fcth UT is given by 

5' + S" 


^IQU° 


\2 • 


(95) 


TV^CX) at (1 + 5)^ 

Furthermore, eonsidering the definition of and in Seetion III-Ai 1, and 

performing simple algebraie operations yields + = 1 + i^n ~ 1) ^in^ p' = J2n=i + 
(e^ — l) sin^ /N, +" = ^ +^n) 1^- This eompletes the proof. 

Appendix D - Proof of Corollary 3 

Applying the matrix inversion lemma [22], and Lemmas 1, 2, and 3, the useful signal power 
of the fcth UT for the IQU-MMSE deteetor in a single-eell uplink massive MIMO system with 
i.i.d. ehannels, IQI only at the UT, and perfeet CSI (i.e., with ^ ^ = iA^i,i,k) is obtained as 


= Jim Pul 




Pul 


AuH 

N 


Ki.kilKlUKi^k 


Pul 


= lim 

(1 + (1 + <5) 


2 ’ 


where A^ = M,k = A^ - Jg.. and 6 is defined as 

ArdAhi.u = %Ltr (D , 


N^qg N 


(96) 


( 97 ) 
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and r is given by 


>iar, 1 v’l 

1 + i ^ IVftjL ) 


Applying similar techniques as in the proof of Corollary 1, we obtain 


5 —^ 


/3^-0 


> -^PulIat- 


(98) 

(99) 

( 100 ) 


Furthermore, the interference power of the /cth UT in the considered single-cell uplink massive 
MIMO system can be expressed as 


rIQU° _ 
^i,k 


K 

^UL 

q=l^q^k 




K 




q=l 


Using the matrix inversion lemma Il22]| and Lemmas [Hm and[3l the first term on the right hand 
side of (IIOII) can be expressed as 

2^1 


q=l^q^k 


q=l^q^k 


a.s. ^ l^pVL^' 


4- 


AT-j-oo (1 + ^) 

Here, 5' is given by 


( 102 ) 


N 




itrfkAlV^ 




\U\"n^pI^. 


(103) 


where we used Lemma 4 and (llOOl) . Following a similar procedure, and considering 5 S> 1, the 
second term on the right hand side of (IIOII) can be approximated as 


K 

Pul 

q=l 




K 


^UKi,kKi,kCA^i 


Pul 




q=l^q^k 


a.s. ^ PUl5^ 




(104) 


N^oo (1 + SY Ca 

Moreover, considering (fT4l) for L = 1, perfect CSI, and no IQI at the BS and applying the matrix 
inversion lemma [|22]| and Lemmas [H [2l and [U the received noise power of the kth UT is given 
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by 


7lQU° _ 


N 






N^oc N{1 + 6YN 


tr (A-^) 


pIl 


TV^cx, iV(l + 5)^ N\U\ 


2 • 


(105) 


Considering (1^ (IIOII) . (11021) (11041) . (11051) . and performing straightforward algebraic simplifica¬ 
tions yields (ISTl) . This completes the proof. 


Appendix E - Proof of Theorem 2 


According to (1421) . for the IQA-WLMMSE detector, the useful signal power corresponding to 
the kth UT in the fih cell is given by 


SlT =0-VuL 


Si i kl 


2N 

2\ 


G.^GL ^ 1 


-1 


-I 


2N 2NprjL 


2N I SiA,k 


-f 


Si,i,k+Kl 


2N 


GmGm + 1 


2N 2iVpuL 


I 


2N 


^ Si,i,k+K 


(106) 


where the first and second terms on the right hand side represent the power of the real and 
imaginary parts of the useful signal, respectively. Applying the matrix inversion lemma l[22]l . we 
obtain the following expression for the first term on the right hand side of (11061) 

I A O' 

2N^i,i,k'^H,kbi,i,k 




2N \ 2N 
where A* ^ is given by 


G^M^ + 1 


I 


2ApuL 


2N 


-1 

2 

\ ~ 


1 Si,i,k 



^ 2N^i,i,k'^i,kSi,i,. 


- 1 : 


A = ^ p p"*" _ ^ ^ 

^,k 2ArpuL 


2N' 


Considering (|3^ and applying Eemmas \T\ [2l and [3] yields 

- 1 , 2 


2N^i,i,k'^i,kSi,i,k 


^ 2N^i,i,k'^i,kSi,i,k 


- 1 ; 


-A 




i.k 


N^oo 


1 + 


2 ’ 


where Xi.fc and ^ are defined as 

, A f 1,0) 


2N 

1 




.(2) 


1 -T • -1 


= lim - 

N^oo 2N 


N^oo 2N 

1 -T ~ 

Sz,Z,/C'^^Z,Z,/cS2,2,/c * 


(107) 


(108) 


(109) 


( 110 ) 

( 111 ) 
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Moreover, $and $are given by 

I O-sE ([“‘.‘Hi + I=«ll) 


i,i,k 


1=1 






and f j is determined as 


1 

m 


K / i(l) 


k=l 




i,i,k ^ ^,^,/c+il 

1 + Si^k 1 + ^i,k+K } 2A^PuL 


+ 


I 


-2N 


where ^f-k+K defined as 

/ L 


^i,i,k+K ~ ^i,i,k | 0-5 ^ ^ [“^.^12,2 

V i=i 


and 5i^k and Si^k+K are the solutions to the following fixed-point equations 

' ^(1) i(i) \ 


A ^ tr I F'‘> f ^ V f , *2,.^ 


2WpuL 


t2Ar 


1 I -,(1) ( 1 f _|_ ^i,i,g+K _ 

^ ' 1 + 1 + ^i,q+K J 2iVP ul 


'5i,k+K I [ 2jy 


■l2Ar 


Similarly, the seeond term on the right hand side of (11061) ean be expressed as 

1 / 1 . . 1 \ -1 ^ 1 ~T A 


2jyk,i,k+K ^ 2NpuL^^^) 


( 112 ) 

(113) 

(114) 


j ^L,k, (115) 


(116) 


(117) 


a.s. 


\Xi,k+K\ 


IN ^i,i,k+K'^i,i,k+K^i,i,k+K 

^ 3~ 2N ^i,i,k+K-^i,i,k+K^i,i,k+K 


N —>-00 /i I ? ^ 

1 J- -h Oi^k+K j 


(118) 


where 5i^k+K is given in (11171) and Xi,fc+x is defined as 

Xi,k+K = {^i,i,k+K^i 


with 


^?lk+K = 0-5 {[^lAlp + i-iAlp) ^^,^,k^iRi,i,k^'^■ 


(119) 


( 120 ) 
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Considering (11061) . (11091) . and (11181) the asymptotie useful signal power can be expressed as 

/ , ,o o \ 


= lim = O.SpuL 


^IQA 


V 


\Xi,kf \Xi,k+Kf 

2 + 


( 121 ) 


1 + Si^kj (^1 + Si^k+Kj j 

According to (321), the interference power of the signal corresponding to the /cth UT observed 


at the output of the IQA-WLMMSE detector of the Ah BS is given by 

L K / „ 2 


^IQA _ 0-5pUL 
“ 4Ar2 




1=1 q=l 
{l,q)^{i,k) 


+ 


- 1 . 


Si,i,k^i Si 




+ 


- 1 , 


Si,i,k+K'^i Si,l,q 


+ 


- 1 , 


Si,l,q-\-K 


( 122 ) 


where Aj = ^Gi^iGj^ + The first term on the right hand side of (11221) can be 


2N 

reformulated as 


L K 

EE 

i=l q=l 
{l,q)jt{i,k) 


1 


E 

1=1 


1 £.J 


L K 

EE 

1=1 q=l 
q^k 


1 ~T A~^~ 

2^j\^Si,i,k-^i Si. 


l,q 


Applying the matrix inversion lemma [l22ll . Lemmas 1-3, and considering (II101) yields 


/ _ ^2) 


\,i,k 


1 + m^,i,kKiki,k 2iv (i + 5.,,) 1 + kk ’ 


(123) 


(124) 


where Xi^i^k is defined as 


Xi,i,k - — tr , 


-( 2 ) 


(125) 


- ^2) 

and 1 ^ is given by (II131) . On the other hand, for the second term on the right hand side of 
(I123L using matrix the inversion lemma Il22ll and Lemmas 1-3 yields IfT^ 


J_-T X-l~ 

2jY 


^ Q*"^ A o’ 

2N ^i,i,k i,k 


1 + mU^,k\kS 


-1; 


i.i.k 


~T *-1^(1) i-i- 


4Ar2 1 + Si^k 


1 + Si j 


I i A.T A“^ ^0) ~ 

^ \ 4iV2 i,k,qSi,l,q + 


2]\l Si^l^q-^i^k,qSi,i,q Si^i,q-^i,k,q^i,i,k'^i,k,qSi,i, 


-1 


-1 iO) I-I 


4Ar2 1 + 


-23fJ< 


( A-P-''’ A ^ \ (—EffT A ^ A ^ 

\2N^i,hq i,k,qoi,l,qj I 47v2oi,Z,q i,fc,g 

1 + 5i,c 


N^oo 


1 + 5^ 


i,k 


0 , 
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2 

S^i,l,q 



1 + 


23fJ< 


1 + ^1 


Qi,l,q 


i.k 


1 + Si I 


2 ’ 


(126) 


where Ai^k,q = JJ.obtained 

by replaeing k with q in (11161) and (11251) . respeetively. Moreover, Ri^i^q, and ^i^i^q are defined 
as 

^ ~T A 5.(^) A ~ a-s- J- ± {a 


~T X-i ^(1) X-i ~ 

Si,l,g i,k,q i,i,k i,k,qSi,l,q ^^2 


Am 

1 CT X-1 ^(1) X-1 

Am^ 


> ^tr f R,,.„f' , ) = Ci,;,, 


Si^i,q-^i,k,q^ i,i,k-^i,k,qSi,i,q 


a-s. 1 , / i (1) fS' 

TTT;Ttr , uT 


i,i,k^ i,k I — '^i,^,q 


^ N^oo Am 

Si,l,q-^i,k,q^i,i,k-^i,k,qSi,i,q ^ ~a1\T2^^ ( ) 4^i,i,q) 


1 ~T I-l 5.(1) I-l 


(127) 

(128) 
(129) 


where f ^ ^ is given by T' in Lemma |4] with N being replaeed by 2N and T, C, and A^. 
being equal to f*, and $-^j^^/(2A^), respeetively. Performing a similar proeedure for the 


other terms in (11231) . the interferenee power as given in (1481) is obtained. Next, we derive the 
asymptotie value of the noise power for the IQA-WLMMSE deteetor. Aeeording to (l42l) . the 
noise power eorresponding to the fcth UT in the ith eell at the output of the IQA-WLMMSE 
deteetor is given by 2 

1 i J A ~1 T ~ 

-^giik-^i ^i'^i 


y^IQA _ 

^i,k — 


2 AT' 


+ 


2N 


^i,i,k+K'^i 


A. ^^,z. 


(130) 


Applying the 


-1 


: T T -1' 


; 7 lQA _ A,fcgi,i,fc ^■S^^i^k+K^i,k+K^i^i^i,k+K^i,i,k+K 

} m TETT ^2 7 “ m I TT- 


-1 


(131) 


2N^i,i,k'^i,kSi,i,kj i“ 2N^i,i,k+K'^i,k+KSi,i,k+K 

Einally, applying Eemmas \T\ [H and [3] to (11311) leads to (l49l) . This eompletes the proof. 


Appendix E - Proof of Corollary 4 

Applying the matrix inversion lemma ll22ll . and Eemmas 1-3 in (11061) . the asymptotie power 
of the useful signal of the kth UT in a single-eell uplink massive MIMO system with perfeet 


CSI and IQI only at the BS is given by 


qIQA°_ 
‘^i,k ~ 


O.SpuL 


- 1 , 




O.SpuL 


- 1 , 


+ 


-1 


Si,i,k+K-^i,k-\-KS\ 


i,i,k-\-K 


Am [1 + A,fc ) 4 Ar 2 M + jj^gi^i^k+K\k+KSi, 


-1 


ii,i,k+K 


PvlS^ 

1 + 5 


2 ’ 


( 132 ) 
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where 




3=—tr ($+f), 

(133) 

and r is given by 

l + ^+2/V^7“j ■ 

(134) 


Assuming /3 <C 1, 5 S> 1, and applying Taylor series expansion, (11341) can be rewritten as 

(135) 

1 + C) 


Substituting (11331) into (11351) and considering /i <C 1, the following quadratic equation is obtained 


+ (1 - 2A^pul/w) S + (4A^VuL/^At - 2Npui^fl) = 0, 


(136) 


straightforward simplifications leads to 


where p = tr [ ) / {2N) and p = ti [ ( ) ) /(2A^). Solving (11361) and applying further 

> 2Npui^p. 


S = 


p+ \^p‘^ + A/3p 


2 V y ^ 

Substituting (11371) into (11351) and considering /S <C 1 yields 

f = 2A^puLl2Ar- 


(137) 


(138) 


Now, we evaluate the interference power of the fcth UT in the considered single-cell uplink 
[MO s> 

O.SpuL 


massive MIMO system, which can be expressed as 

2K 


rIQA°_ 

'^i,k 




2K 


Am 


~T A ~ 


+ 


- 1 , 




1,1,q 


2K 


Pul ~t a“1~ 

3jY2 i Si,i,qSi,i,q i Si,i,, 

q=l,q^k 


Pul 

8Ar2 


^li,k+KK ^Zi,i,qZli,qK 


2,2,fc+i^ 


q=l,q^k 


a.s. ^ PVlI^S' 
N^oo 


1 + 5 


4 ’ 


(139) 


where we applied the matrix inversion lemma Il22l and Lemmas 1 and 2. Here, S' is given by 


S' = ^tr ( 4'#'^A, ^4'4'^A, ^ 


^tr i t 


(ms) 


^ -^tr ( 


(140) 


where we used Lemma 3. Assuming (S <^1 leads to 

S' = 4Ar2p2 


( 141 ) 
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Moreover, the reeeived noise power eorresponding to the fcth UT is given by 


7lQA° 

^i,k 


-1 




+ 


-1 






1 a~^t~~Tx'*'a~^~ 

Si,i,k+K^i A. g 


4Ar2 


-A 


4Ar2 

S' 


i oi i k+K ~ / \ 2 ’ 

Arfl + 5 


(142) 


where we applied the matrix inversion lemma [l22ll and Lemmas 1-4. Now, eombining (11321) . 
(11391) . (11421) . and 5 S> 1 yields 


SINK: 


IQA° 

i.k 
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Next, we derive elosed-form expressions for p = tr /2N and p = tr j /2N. 

It ean be easily shown that nil^ = l 2 Ar. Considering this and using ([3]), we obtain p = 


En=il|A„||y2iV = 1 + Similarly, p is obtained as p = ^ En=i (1 + + 


1 ^2 
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(1 — en)^ -I- 2 (1 — e^)^ sin^ 6'„. Using e„ <C 1 and performing simple algebraie operations yields 
p = 1 -f ^ (12 — 4sin^ 6n) -f 2 sin^ 6*„. This eompletes the proof. 


Appendix G - Proof of Corollary 5 

From (11061) . the matrix inversion lemma ll22ll . and Lemmas 1-3, we obtain the asymptotie 
power of the useful signal of the fcth UT in a single-eell uplink massive MIMO system with 


perfeet CSI and IQI only at the UT as 
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where 


~ + [“Lfc]2,l f-fP\ 
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Applying similar teehniques as in the proof of Corollary 4, we obtain 
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Substituting (11471) and (11481) into (11441) yields 
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Performing straightforward algebraie manipulations leads to 
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Applying similar teehniques as in the proof of Corollary 4, the asymptotie interferenee power 

ystem ean be expressed as 

_j ~T 


of the Pth UT in the eonsidered single-eell uplink massive MIMO system ean be expressed as 
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Similar to the proof of Corollary 4, the asymptotie noise power of the Pth UT ean be written 
as 

1 -X T-1. 
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Observing (11521) and (11531) and considering /? <C 1 and 5k, 5k+K 3> 1 the interference power 
vanishes and the only remaining disturbance is the noise power. Taking this into account and 
performing straightforward algebraic operations leads to (|58]) . This completes the proof. 
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